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^ ■ Abstract 
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An intrinsic multitype branching structure within the transient (1, i?)-RWRE 
is revealed. The branching structure enables us to specify the density of the 
absolutely continuous invariant measure for the environments seen from the par- 
ticle and reprove the LLN with an drift explicitly in terms of the environment, 



comparing with the results in Bremont (2002). 
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§ ! 1 Introduction and main results 



A random walk in random environment (RWRE, for short) {X n ,n 6 N} with bounded jumps 
on the line, written as (L, 72)-RWRE, means that for each step the possible jump range to the 
left is bounded by L and to the right by R, where L and R are positive integers. The aim of this 
paper is to reveal the branching structure within the (l,i?)-RWRE, which is a story different 
with but essentially complement for that of (L, 1)-RWRE (Hong and Wang, [10], 2009) when it 
transient to the right. 

It is well-known that when the walk is transient, i.e., X n — > oo, the intrinsic branching 
structure within the (1, 1)-RWRE is a Galton- Watson branching process with geometric offspring 
distribution, which plays an important role in the proofs of the limiting stable law (Kesten et 
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al [12], 1975), the renewal theorem (Kesten [llj . 1977), and the law of large numbers (Alili pQ, 
1999) for this nearest random walk in random environment, see also Gantert and Shi ([9], 2002). 

However, when L or R > 1, and the random walk is transient (we consider X n — > oo), 
It seems no result about the intrinsic branching structure within the (L, i?)-RWRE up to our 
knowledge, even for random walks in non-random environment. Recently, partially progress has 
been made by Hong and Wang ([10] . 2009), where a multitype branching structure within the 
(L, 1)-RWRE have been revealed when X n — > oo and the Kesten's type stable law have been 
proved. In the present paper, we focus on the opposite direction which is more complicated: we 
will figure out the intrinsic branching structure within the (1,-R)-RWRE when X n — > oo(note 
that the situation for the (L, 1)-RWRE transient to the left is equivalent to that for the (1,12)- 
RWRE transient to the right). We will discuss R = 2 in detail and extend to general R at the 
end. 

RWRE has been studied extensively in recent years, especially for the nearest (1, 1)-RWRE 
on which many results have been obtained, see for example [18] . |12j . |llj . [2], etc., we refer to 
Sznitman ([IT], 2002) and Zeitouni ([J9], 2004) as a general review. For the (L, J R)-RWRE, Key 
(|13j. 1984) discussed recurrence/transience criterion in terms of the sign of two intermediate 
Liapounov exponents of a random matrix. Bremont ([6], 2009) formulate a criterion for the 
existence of the absolutely continuous invariant measure for the environments seen from the 
particle and deduce a characterization of the non-zero-speed regime of the model. It should 
note that the (L, 12)-RWRE can be treated as a special case of the random walks in random 
environments on a strip, see Bolthausen and Goldsheid ([3], 2000; [I], 2008 ) and Roitershtein 
([H], 2008). 

Bremont ([5], 2002) has proved a recurrence/transience criterion for (L, 1)-RWRE involv- 
ing the greatest Lyapunov exponent with respect to a random matrix M , and the law of large 
numbers by assuming the (1-M) condition related the existence of an invariant measure for the 
environments seen from the particle. As an important application of our branching structure 
within the (1, 1?)-RWRE, it enables us to specify the density of the absolutely continuous invari- 
ant measure explicitly and reprove the LLN with an drift explicitly in terms of the environment. 

We now introduce the model of (1, 2)-RWRE. Generally speaking, random walks in random 
environment involves two kinds of randomness: the transition probability (we call it the "envi- 
ronment"), which is chosen from a specified distribution; and the random walk driven by the 
chosen "environment". Specifically, let A = { — 1, 1, 2} be the set of possible jump range of the 
random walks. Let A^(A) be the collection of all probability measures on A. Then define an 
environment to be an element oj = {(q(oj) x ,pi(oj) x ,p2(oj) x ) : x G Z} G A / f(A) z =: O. Let P 
be a stationary and ergodic probability distribution on (£l,P) and 9 be the spatial shift, i.e., 
(0ui) n = uj n +i- Assume that 9 is an invertible transformation on the probability space (f2, P), 
measurable as well as its inverse and preserving P. Moreover, assume that the environment is 
elliptic: 

3e>0, V z G {1, 2}, (p z /q) > e, P-a.s.. 

Given an environment uj G Q, one can define a random walk {A n } in the environment oj to be 
a time-homogeneous Markov chain on Z with Xq = and the transition probabilities 

P w (A n+ i = x - l\X n = x) := q(9 x u) = q(u) x , 
P w (A n+ i = x + z\X n = x) := p z (9 x uj) = p z (uj) x , 
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for all x £ 7L and z £ {1, 2}. For each cj, we use P u to denote the law induced on the space of 
paths (Z N , Q). Then, define a probability measure P := P <g> P w on (0 x Z N , J" x 0) by 



P(F x G) := y P w (G)P(dw), F £ J 7 , G £ Q. 



Statements involving P w and P are called quenched and annealed, respectively. Generally, with 
a slight abuse of notation, P can also be used to denote the marginal on Z N . Expectations under 
P w and P will be denoted by and E, respectively. 

We will need the following notations. p z (0 x uj)o will be simply denoted by p z (x). Any 
expression of the form f(9 k co) will be simply denoted by f(k). For the random walk {X n } and 
k £ Z\{0}, we write for the quenched probability starting at k, and E% for the corresponding 
expectation. 

Assume that Xq = and the (l,2)-random walk is transient to the right, i.e., X n —> oo, 
P-a.s.. Let To = 0, and 

T k =M{n>T k _ 1 :X n >X Th _ 1 }, k > 1 

be the sequence of ladder times of the random walks. Note that T k < oo, P-a.s.. To calculate 
T\ accurately, we will figure out a multitype branching processes by decomposing the path of 
the walk. Intuitively, if the walk from i < take a step to i — 1, it must crossing back to i 
or jumping over i (to i + 1) because of T k < oo, P-a.s., in which there are only three kind of 
backing ways: from i — 1 to i, from i — 2 to i and from i — 1 to i + 1. So we divide all the steps 
from i to i — 1 into three kind of steps according the crossing back ways. Let A(i),B(i), and 
C(i) are the numbers of steps from i to i — 1 before time Ti with crossing-back from i — 1 to i, 
i — 2 to i and z — 1 to i + 1, respectively. And for the last step of Ti, we can consider it as a 
immigration for the multitype branching processes. 
Set for i < 0, 

U(i) = [A(i),B(i),C(i)}, 
Then we have the branching structure within the (1, 2)-RWRE as follows. 



Theorem 1.1 Assume X n oo, P-a.s.. Then for P-a.s. uj, (U(i) = [A(i),B(i),C(i)]j is 
an inhomogeneous multitype branching process with immigration 

U{1) = [1,0,0], with probability 

U(l) = [0,1,0], with probability 

U(l) = [0,0, 1], with probability 

The offspring distribution is given by 

P u (ll(i) = [a,b,Q] | U(i + 1) = [1,0,0]) = [1 - a(i) - f3(i)]C^ +b a(i) a (3(i) b , (1-1) 
P u (U(i) = [a, b, 1] | U(i + 1) = [0, 1, 0]) = [1 - a(i) - /3(z)]C a a +fc a(i) a /3(i) b , (1.2) 



Pi(0) 


1 - a(0) - 


/3(0)' 


7 (0) 




1 - q(0) - 


mv 


p 2 (o) 




1 - a(0) - 


P(0) 
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P u (U(i) = [a, b, 0] | U(i + 1) = [0, 0, 1]) = [1 - a(i) - P(i)]C a a+b a(i) a (3(i) 1 ', (1.3) 



where 



1 (i)=q(i)-P^ 1 [(-oo,i-l),i + l], 



a(i) =g(i)-Pr i [(-oo,i-l),i] 

/3(i) = (7 (i)-pr i [(-oo,i-i),i] 



- 1) + 7(i - 1) ' 
pi(i - 1) +7(i - 1)' 



and P^ oo, i — 1), z — 1 + j] ■ = P^ 1 {reach [i, +oo) for the first time at the point i — 1 + j} 
for j = 1,2, the exit probabilities which can be expressed in terms of the environment in Lemma 
IP □ 



As an immediate consequence of Theorem 11.11 We can get the offspring quenched mean 
matrix of the multitype branching process and the quenched mean of the T\ . 

Corollary 1.1 Assume X n — > oo, F-a.s.. Then for P-a.s. u, and i < 0, the offspring mean 
matrix of the {—i + l)-th generation of the multitype branching process is 



( 



N(i) 

Moreover, 

E U (T{) = 1 + ((2,2,1), 



l-a(i)-f3{i) l-o(i)-/3(j) 



i 



l-a(i)-/3(i) l-a(i)-/3(i) 
\ l-<*(i)-/3(i) l-a(i)-P(i) U / 



1 



1 - a(0) - /?(0) 



P1 (o), 7 (o), pa(o)) • J3j\t(o)--.jv(o). 



j<0 



□ 



As an application of the branching structure to the random walks in non-random environ- 
ment, let 

X = 0, X n = £i H + £n, 

where £i, £2, • • • is a series of i.i.d random variables with 

P(6 = -l) = ?, P(6 = i)=pi, P(6 = 2)= P2 . 

The computability of E(T\) by the branching structure as in Corollary 11.11 enable us to validate 
the Wald's equality. 

Proposition 1.1 Assume that P(£i) = Pi + 2p 2 — q > 0. Then the Wald's equality holds: 

E(X Tl ) = E{T X ) ■ E{X X ). 



□ 
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From the point of view "environment viewed from the particles", define uJ(re) = 6 x "uj. Then 
the process {U(n)} is a Markov process under either P w or P, with state space f2 and transition 
kernel 



e 2 a 



Set ip\ kw = P ekuJ (X Tl = 1), and tp* hu = P ekul {X Tl = 2) = l-y\ hu . Whenever E(Ti) < oo, define 
the measure 

lx^=2^\ \ Q(5) 



/l -1 Tl_1 1 -2 Tl_1 \ 

Q(5) = El ^ l{s(i)es} + ^We-B} I. <2( 5 ) 



Q(«) 



The significance of the branching structure to the (1, 2)-RWRE is that we can express the density 
of Q with respect to P explicitly. 

Theorem 1.2 Assume E(Ti) < oo. Then Q(-) is invariant under the Markov kernel K , that 
is, 



Q(B) = J J l^ B K{u},du/)Q{du). 



Furthermore, 



where 



dQ 
dP 



n(w), 



2 + ((l,l,l),Ei 



Pi(») 7(0 



*>1 Vpi(i)+7W ' PiW+TW ' 



1 ■ N(i) ■ ■ ■ N(l) 



1 - a(0) - /3(0) 



□ 



Therefore, the LLN for the (1, 2)-RWRE can be reproved with an explicit drift. 
Theorem 1.3 Assume E(Ti) < oo. Then 

lim — - = vp, P-a.s., 

n— >oo TL 



where 



Ei 



Vp 



Pi(0)+2p 2 (0)-P-i(0) 
X-ol{Q)-P(0) 



2 + ( (1, 1, 1), E,>i ( pTiSw ' pTlJrW 1 ) • * « • • • 



2 + ((2, 2, 1), E t > i^ify ^Swy l) ■ * (0 ' ■ ■ * (0) 



□ 
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We arrange the remainder of this paper as follows. In section [2j we formulate the intrinsic 
multitype branching structure within the (1,2)-RWRE under the assumption X n — > oo, i.e., 
Theorem ll.lt Corollary 11.11 and Proposition 11.11 are proved. In section [3l we specify the density 
of the absolutely continuous invariant measure for the environments seen from the particle 
explicitly based on the branching structure, i.e., Theorem 11.21 is proved. In section [H as an 
application of the branching structure, the LLN will be reproved with an explicit drift, that is, 
Theorem 11.31 is proved. Finally in section [5j for general R> 1, we give the intrinsic multitype 
branching structure within the (1, i?)-RWRE under the assumption X n — > oo. 



2 Intrinsic multitype branching structure within the transient 
(1,2)-RWRE 

We introduce the following exit probabilities of leaving a given interval from the right side. 
Consider integers a,b,k with a < k < b, define 

P^[(a, b), b + 1] = P^jreach (b, +oo) before (— oo, a) and at the point b + 1}, 
P^[(a, b),b + 2} = P^{reach (b, +oo) before (-oo, a) and at the point b + 2}. 

These exit probabilities play an important role in the offspring distribution of the branching 
structure, which can be expressed in terms of the environment (Bremont [5J, page. 1273-4, Lemma 
2.1 and Proposition 2.2). Only a slight modification should be made: exiting from the left in [5] 
corresponds exiting from the right here. We still give the details of the proof for convenience. 

Lemma 2.1 For n > 2, we have 

( ei , [M(i) + ■■■ + Af(-n) • • • M(i)]v) 



I»[(-n,i),i + l]: 
P*[(-n,i),i + 2] 



1 + {e x , [M(i) + • • • + M(-n) • • • M(i)] ei ) 
(ex, [M(i) + ■■■ + M(-n) ■ ■ ■ M(i)]e 2 ) 



1 + (ex, [M(i) + ■■■ + Af(-n) • • • M{i)]ex) ' 
where ex = (1, 0)', e 2 = (0, 1)', v = ex — e 2 and 

( Pi(i)+P2(i) P2(i) \ 

Furthermore, if X n — > +oo, P-a.s., then 

P^[(-oo,i),i + l] + P^[(-oo,i),i + 2] = l. (2.1) 

Proof. Set f(k) = P^{(—n,i),i + 1]. For k such that — n < k < i, using Markov property, we 
get a harmonic-type recurrence equation: 

f(k) = px(k)f(k + 1) + P2 (k)f(k + 2) + q(k)f(k - 1), 



G 



with f(i + 1) = 1, f(k) = for k > i + 2 or k < -n - 1. Setting g(k) = f(k) - f(k - 1), we 
obtain 

9{k) = ^ 9{k+i)+E M( 9{k+2)+9{k+i) 



q[k) q(k) 



Rewrite the equation as 



g{k) x = / ejfi\, g(k + 1) 

g(k + i) J { l o J\g(k + 2) 

Set 

^ - ( ) ■ ^ ^ - ( °f > °f ' 

where ai(fc) = Pl ^ k g^ 2 ^ , a>2(k) = p g 2 / fc 7 ■ Thus, we obtain the relations 

[/■(fc) = M(k)U(k + 1) 

= Af(ife)---M(i)t7(i + 1), 

with C/(i + 1) = ([1 - /(i)], -1)' and U(-n) = (/(-n), [/(-n + 1) - /(-n)])'. Summing from 
— n to i, we deduce that 

/(») = (e l ,M{i)U{i + l)^ + --- + (e l ,M{-n)M{-n + l)---M(i)U{i + l)^ 
M(i) + ■ ■ ■ + M(-n) • • • M(i) 



1-/(0 
-1 



The first formula then follows. By similar reasoning, one can get the second formula. The second 
part of the lemma is just the conclusion of Proposition 2.2 in [5]. □ 

Now we introduce the branching structure. Recall Xq = and the sequence of ladder times 
of the random walks: Tq = and 

T k = inf{n > T fc _! : X n > X n _J, k > 1. 

Note that T k < oo, P-a.s. if X n — > +oo, P-a.s.. Define, for i < 0, 

T7i )0 = min{£; < Ti : X k = i}, 

9 i>0 = min{77i i0 < k < T x : X k ^ = i,X k =i-l}, 

and for j > 1, 

Oij = mxn{9ij-i < k < T x : X fc „i = i - 1, X k = i}, 
(3 it j = min{6'jj_i < k < T\ : X fc _i = i - 2,X fc = i}, 
7ij = min{6'ij_i <k <T\\ X k _ x = i - l,X k = i + 1}, 
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9ij = mm{r]ij < k < T x : X fc „i = i,X k = % - 1}, 

with the usual convention that the minimum over an empty set is +00. We refer to the time 
interval as the j-th excursion from i—1 to {i, i + 1}. For any j > 0, any i < 0, define 

A,j = #{k > : Oi+ij < 9 i)k < r] i+1 j +1 , and r/ ijfc +i = "i.fe+i}, 
Bi,j = #{/ > : i+ ij < 0^ < Tft+ij+i, and 77,^+1 = A,z+i}, 
Cij = #{m > : 9 i+ ij < 9 i)m < r] i+1)j+ i, and r/i >m+ i = 7i, m +i}- 

Note that Aij,Bij, and Cjj are the numbers of steps from i to i — 1 during the (j + l)-th 
excursion from i to {i + 1, z + 2} with crossing-back from i — 1 to z, i — 2 to i and i — 1 to z + 1, 
respectively. Special attentions should be paid to the different ending ways of each excursion 
and their consequences, take the (j + l)-th excursion from i to {i + 1, i + 2} for example: 

• If T]i + ij + i = CKj+ij+i, then the excursion [0j +lj - +1 , aj +1 j +1 ] ends by jumping from i to 
i + 1. " 

• If T/j + i J+ i = 7j + ij + i, then the excursion [0j + ij + i, 7j + ij + i] ends by jumping from i to 
i + 2. 

Furthermore, in the above two cases, excursions from i — 1 to {i, i + 1} included in the (j + l)-th 
excursion from i to {% + 1, i + 2} must end at i and consequently, = 0. 

• If r]i + ij + i = /3j + i J+ i, then the excursion [0j + ij + i, /3j+ij + i] ends by jumping from i — 1 to 
i + 1. ' 

Moreover, in this case, the last excursion from i — 1 to {i, i + 1} included in the j + 1-th excursion 
from i to {i + 1, i + 2} also ends by jumping from i — 1 to i + 1, that is, if 9i )TnQ is the beginning 
of it, then 77^+1 = 7i, mo +i = ft+ij+i = r/i+ij+l- Therefore, Cjj = 1. 

Define for i < 0, 

aw = 2 ^, = E B ^ c « = E c m 

j>0 J>0 j>0 

to be the numbers of steps from i to i — 1 before time Ti with crossing-back from i — 1 to i, i — 2 
to i and i — 1 to i + 1, respectively. Define for i < 0, 

l7(i) = H(i),£(i),C(i)]. 

Then A{i) + -B(i) + C(i) is the total number of steps the walk jumping from i to i — 1 before 
T\. The number of steps crossing back from i — 1 to i is A(i) + B(i), because X n — > 00, P-a.s. 
and C(i) = B(i + 1) is counted as steps crossing back from i to i + 1. Thus we have 

Ti = l + ^(24(t)+2B(t) + C(i)) 

i<0 

= 1 + ((2, 2,1), $>(<))■ (2-2) 

i<0 
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In order to study Ti, we consider {£7(i)}j<o instead. We need the following probabilities. For 

i < 0, 

a(i) = P^[ the walk jumping from i to i — 1, and finally jumping back (2-3) 

from i — 1 to i before Ti], 

(3(i) = P^[ the walk jumping from i to i — 1, and finally jumping back (2-4) 

from i — 2 to % before Ti], 

= P^[ the walk jumping from i to i — 1, and finally jumping back (2.5) 

from i — 1 to i + 1 before Ti]. 

Expressions of these probabilities will be given later in terms of the exit probabilities. By Markov 
property, and the different ending cases explained above, we obtain, for % < 0, 

P u (u(i) = [a, 6,0], U(i + l) = [1,0,0] | |l7(t + l)| = l) = CZ +b a(i) a /3(i) b Pl (i), (2.6) 

P w (?7(i) = [o,6 > l],t7(i + l) = [0,l > 0] | |l7(t + l)| = l) =C a a+b a(i) a p(i) b 1 (i), (2.7) 

P w (l7(i) = [a,6 ) 0],^(i + l) = [0,0,l] | |^(i + l)| = l) = C a a+b a{i) a ^{i) b P2^ ■ (2.8) 

Additionally, for i = 0, we define £7(1) = [^4(1), B(l), C(l)] for consistency. Noting that the walk 
starts at 0, and the excursion from to {1,2} ends at time T\, there is only one jump crossing 
up from to 1 in the time interval [0, T\], which can be regarded as |J7(1)| = 1 if T\ < oo, P-a.s.. 
If the excursion from to {1, 2} ends by jumping from to 1, then set A(l) = 1, in other words, 
U(l) = [1,0,0]. By similar reasoning, £7(1) = [0,1,0] if the ending jump is from —1 to 1 and 
£7(1) = [0, 0, 1] if the ending jump is from to 2. Summing over a, b > in (12 .6p . we have 

P u (u(i + 1) = [1,0,0] I |l7(t + l)| = l) 
= E P «( U ® = [a,b,0],U(i + l) = [1,0,0] | \U(i + l)\ = 1 



a,b>0 



2 CZ +b a(i) a m b Pi(i) 



a,b>0 



l-a(i) -/9(t)' 



(2.J 



Similarly, 



P u (U(i + l) = [0,1,0] |£7(i + 1)1 = 1 
P w fl7(t + 1) = [0,0,1] |£7(i + 1)| = 1 



7(i) 



(2.10) 
(2.11) 



1 - a(i) -/3(i)' 

P2(») 

1 - a(i) - ' 

Now we are ready to calculate a(i), f3(i), and 7(2). Firstly, by the definitions (12.3p - (l2.5jl 
we can see 



a(i) = P u (u(i) = [1,0,0] I \U(i)\ = l), 
P{i)=P u {u{%) = [0,1,0] I |C/(i)| = l), 
7(i) = P w («7(i) = [0,0,1] I |£7(i)| = l). 
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Thus from (12. 9p and (|2.10p . we have the ratio 

a(i) : = pi(i - 1) : j(i - 1). 
On the other hand, by the definitions (|2.3p - (|2.5p again, we know 



(2.12) 



a(i) + (3(i) = P^[ the walk jumping from i to i — 1, 

and finally jumping back to i before T\] 
7(i) = P^[ the walk jumping from i to i — 1, 

and finally jumping back to z + 1 before Ti] 

Recall the exit probabilities, we obtain 

a(*)+0(i) = g(i)-Pr 1 [(-oo,i-l),i], 
7 (i) = g(i)-i^ _1 [(-oo.i- + 



Finally, combining f|2. 12|) and (|2.13p we get 

a(i) = (? (i)-P ( :- 1 [(-oo,i-l),z]- 
/3(i) =?(*)■ *-!),*] • 



- 1) 



- 1) + 7(i - 1) ' 
Pi(i - 1) +7(i - 1)' 



(2.13) 
(2.14) 



(2.15) 
(2.16) 



Proof of Theorem[TJ\ By and (l2~^j) - (j2TTTD . we obtain 

P w ([/(f) = [a, 6, 0] | I7(t + 1) = [1, 0, 0]) = [1 - a(i) - p(i)]CZ +b a(i) a {3(i) b , (2.17) 

Pu (U(i) = [a, b, 1] | U(i + 1) = [0, 1, 0]) = [1 - a(i) - /3(i)]C a a +b a(i) a /3(i) 6 , (2.18) 

P w (0 = [a, 0] | + 1) = [0, 0, 1]) = [1 - a(i) - /3(i)]C a \ b a(i) a /3(i) fe . (2.19) 



(|2.17p - (|2.19p give the offspring distribution for the (— i + l)-th generation (i < 0). Indeed, it is 
enough to show 



E 

a,6>0 



PjC/(i) = [a,6,0] tf(i + l) = [l,0,0])-PjE/(i + l) = [1,0,0] |C/(i + l)| = l 



+^(t/(i) = [a, 6,1] I7(t + l) = [0,l,0]j ■P u [U(i + 1) = [0,1,0] |CT(i + l)| = l 
+i' w ( , ?7(i) = [a,6,0] 17(* + 1) = [0,0,1]) • pj[/(i + 1) = [0,0,1] |CT(i + l)| = l 



which is equivalent to check 



Pi(i) + +P2(i) 



1. 



l-a(i) - P(i) 

By (|2.ip . (|2.13p . and (|2.14|) . we know a(i) + f3{€) + = q(i). The conclusion therefore follows 
from pi(i) +P2{i) + q{i) = 1- 
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Since Xq = and X n — > oo, P-a.s., we have T\ < +00, P-a.s., so we imagine that there 
are different types of particles immigrating into the system U(l) = [A(1),B(1), C(l)] and hence 
\U(1)\ = 1 P-a.s.. By (|2J|, we know 

Pi(0) 



P w (tf(l) = [l,0,0]J = P w (tf(l) = [1,0,0] |t/(l)| = l 
By similar argument, 

p„(c/(i) = [0,1,0]) = - 

p.fc/a) = [0,0, 1] 



1 - 0(0) - m • 



7(0) 



«(0)-/3(0)' 

P2(0) 



1 - a(0) - /3(0)' 

Thus the immigration of the multi-type branching process follows. This completes the proof. □ 



Proof of Corollary Summing over b > in (12. 17ft . we obtain 
P w (A(z) = a J7(* + 1) = [1,0, Op 



So the expected number of type-A offspring of a single type-A particle in one generation is 

«(*) 



a>0 



iVn(i) = £VP w L4(i) = a U(i + 1) = [1,0,0] 



l-a(i) — " 



By the same argument, we can get all the elements of the offspring mean matrix N(i). For the 
second result, noting that by the multitype branching process, we have 



and 



E u [U(i)] = E u [U(i + 1)] • N(i) = E U [U(1)] ■ N(0) ■ ■ ■ N(i), 
E U [U(1)] 

= [1,0,0] -P U (U(1) = [1,0,0]) + [0,1,0] -P u {lJ(l) = [0,1,0]) 

+ [0,0, 1] -P u (u(l) = [0,0,1]) 

Pi(0) 7(0) P2 

1 - q(0) - (3(0) ' 1 - a(0) - /3(0) ' 1 - a(0) - (3(0) 

pi(0), 7 (0), p 2 (0)Y 



1 _ Q ( ) _ (3(0) 

The desired conclusion follows by taking the quenched expectation in (|2.2p . 



□ 



As an immediate application of the branching structure, consider random walks in non- 
random environment. Let £i,£2 ; " • " be independent variables with common distribution 



p(& = -!) = <?, p(ei = i)=pi, p(6 = 2)= P2 . 
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The induced random walk is the sequence of random variables 

Xq = 0, X n = £1 H + £n- 

Define Ti = inf{n : X n > 0} as before. The countability of Ti in Corollary 11.11 enable us to 
check the Wald's equality ([8], pg. 397) directly. 

Proof of Proposition First, we have 

E(X Tl ) = 1 • P(X Tl = 1) + 2 • P(X Tl = 2). 

By Lemma 2.1, we obtain 

P(X Tl = 1) = i*[(-oc,0>, 1] = to ^±^±-^ig ,2.20) 

n->oo 1 + (ei, [M + M z + • • • + M n Jei) 

P(X Tl = 2) = P»[(-oc,0),2] = lb, - lei / lM ,-t M lt'' +M 2i , ■ (2.21) 
v 1 ; LV ; J n^oo l + ( ei ,[M + M 2 H hM n ]ei) v ; 



where 



P1+P2 P2 

M := I ^ c ' 
1 



Note that 



1 1 J \ y Ai - A 2 V -1 Ai 
where Ai )2 = Pi+P2=i= V(pi+p 2 ) + 4, ?P2 arg ^ e e ig enV alues of M. Then 

in+l _ /i , -v \ \n+l 



(ei,M"«) = (1 + A 2 )A? +i - (1 + A!)A" +i , 

V 2 ' 



( ei ,M"e 2 ) = (-A 2 )A^ +1 + A 1 A 2 t+ 
(ei,M n ei) = A™ +1 - A™ +1 . 



Since EXi > 0, we have Ai > 1, and A 2 G (-1,0). Hence, by (|2.20p and (|2.2ip . we obtain 

P(Jf Tl = 1) = 1 + A 2 , 
P(X Tl = 2) = -A 2 . 

Therefore, £(X Tl ) = 1 - A 2 . 

The next step is to calculate E(T\), which is done by the branching process (Corollarv ll.ip . 
Since the environment is not random, there is no site index or integration with respect to P any 
more. Applying Corollary 1 1.1 1 to this random walk, we have 

e{t x ) = 1+ ((2,2,1), __L_( P1) 7 , P2 yj2 Nn )- 

a P n>l 
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where 



N :-- 



a 







1-a- 


-0 




a 




l-a- 


-P 


l-a-B 


a 


P 



Using the eigenvalues and eigenvectors of N, and noting that the norm of the greatest eigenvalue 
of ./V is less then 1, we have 

fa p p \ 

H Nn = ^~o ™ \ 2a 2 ? l"2«-/3 • (2.22) 

Thus, 

E(Ti) - 7 + 1-0-/3 



(I- a- /3)(1 -2a- 3/3)' 
Substituting 7 = -q\ 2 , a = q(l + A 2 ) pi ^ A2 , and /3 = q(l + A 2 ) p ~^ 2 , we have 



7 + 1 - a - /3 = A 
1 
2 



l-a-0 = -(!-? + A), 



x _ = (1 + ^-^-3^ + 3^ 

2( Pl -p 2 + A) 



where A = y (pi + p 2 ) 2 + 4p 2 g. Then 



2A 



and 



/-, „, A \ (l-g+A)(l+g-A)( Pl -3p 2 +3A) 
V 1 <? + ^J 2(pi-p 2 +A) 

2A 

f-i _ „ 1 A \ _ 2q(pi -p 2 + A) (pi -3p 2 +3A) 
^ 5 + ^ 2(p!-p 2 +A) 

2A = _ a(p 1 ,p2,q) 

1 - Q ~ QPi + 3<?p 2 + (1 - 3g)A ' b{p 1 ,p 2 ,q) ' 

gtgjQ = 1-A 2 = 3g - 1 + A = _ c{ Pl ,p 2 ,q) 
E(Xi) pi+2p 2 -q 2q(p 1 +2p 2 -q) ' d(p 1 ,p 2 , q) ' 

In order to prove E(Xt 1 ) = E(Ti) ■ E(X\), it suffices to show 

a(pi,P2,q)d(pi,p 2 ,q) = b(p 1 ,p2,q)c(p 1 ,p 2 ,q). 

In fact, 

a(pi,P2, q)d{pi,P2, q) = 4g(pi + 2p 2 - 9) A 

= 4g(l-p 2 -g + 2p 2 -g)A 
= 4g(l +p 2 - 2g)A. 



13 



And, 

b(pi,P2,q)c(px,p2,q) 
= [(1 - q ~ QPi + ZqP2) + (1 - 3?) A] (3? - 1 + A) 

= (1 - q ~ qpi + 3gp 2 )(3g - 1) + (1 - 3<7)A 2 + [(1 - q - q Pl + 3qp 2 ) + (1 - 3g)(3g - 1)]A 
= (3q - 1)(1 - q - qpi + 3qp 2 - (1 - <?) 2 - 4gp 2 ) 

+ [1 - q - q(l - q - p 2 ) + 3qp 2 - (1 - 6g + 9q 2 )]A 
= (3q - l)(q(l -pi- p 2 ) - q 2 ) + [4q + 4qp 2 - 8q 2 ]A 
= 4q(l+p 2 -2q)A. 

Thus, the desired conclusion follows. □ 
Remark Proposition 11.11 is a strong evidence to validate the branching structure. 



3 Density of the absolutely continuous invariant measure 

Now we introduce the machinery of the "environment viewed from the particle". The first 
step consists of introducing an auxiliary Markov chain. Starting from the RWRE {X n }, define 
ZJ(n) = 9 x "uo. The sequence {uJ(n)} is a process with paths in f2 N . This process is in fact a 
Markov process. The proof is the same as Zeitouni ([19], page 204, Lemma 2.1.18), we omit the 
details. 

Lemma 3.1 The process {uJ(ra)} is a Markov process under either P u or P, with state space Q 
and transition kernel 

K(u, dw') = ql^Q-iu +pxl u '=0u +p 2 l ujl= 02 (JJ . 

□ 

The next step is to construct an invariant measure for the transition kernel K. Assume that 
X n -)• oo, P-a.s., implying T n < oo, P-a.s.. Set = P e k w {X Tl = 1), and Lp 2 kuj = P ekui (X Tl = 
2) = 1 — </?gfc w - Whenever E(Ti) < oo, define the measure 



t r -,-^^ \ aim = 

Note that Q is a probability measure. 



Proof of Theorem — invariant measure. We will show 

Q(£) = ^ VgBifCw.dwOOCda;). 

On one hand, 

°°/l \°°/l \ 

Q(fl) = Y, E ( ^T^; Ti > i; S7(i) G fl 1 + E ( 2\ > i; co(i) G 5 ) . 

i=0 V Voj / j=o V 
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On the other hand, 

J I l ul&B K{u,AJ)Q(duj) 

= ^E( Xt \ =1 ; T x > i; 57(i + l) G PJ +^E( X ^ =2 ; Ti > i; W(i + 1) £ P 

„-_n V ¥>u / i=Q 



i=0 



V tu) / ,7=1 w 

+e(" 1x ^ =1 ; Ti < oo; ZJ(Ti) G b) +e( 1Xti 2 =2 ; T x < oo; u(T x ) G £ 



It only needs to show 



E 



lx 



Ti > 0; w(0) £B =E 



Ti < oo; uJ(Ti) 6 5], i = 1,2. 



Indeed, 



R.S. = p p 



1 



Euj(lx Tl =i, w(Ti) G P, Ti < oo 



P w (57(Ti) G P, Ti < oo 



l w G B • Ti < oo 



= P MPw G PJ (since P w (Ti < oo) = 1). 

Similarly, L.S.= P(w G P). Then, by the invariance of the environment, we get L.S.=R.S.. This 
completes the proof. □ 



Proof of Theorem — density. Let / : ft — > R be measurable. Then, 

,Ti-l -, __ Ti-l 



//dQ = E ( r £/(57(t))%-^+£ /(57(f)) 
7 V i=0 ^ w i=0 



e ( £ /(^^ + £ /(^^) 



1 ^T 1 =2 
lx Tl =2 S 



i<0 TUJ i<0 

where Vi = #{k G [0,Ti) : = i}. Using the shift invariance of P, we get 
fdQ = V E P (f(9 l co) \\E w (Vi ■ 1 X =i) + \E u (Vi ■ lx Tl = 2 } 



^Pp(/(u;)[^— P,-^(^ • lx Tl =l) + -^—Eg-iJVi ■ 1 Xt =2) 

Pp(/H^ l-^—Eg-i^Vi ■ lx Tl =i) + -2^— Ee-i^Vi ■ lx Tl=2 )l) 
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Ep(f{u) Y, [Eo-MXtl = 1) + E e - iu {Vi\X Tl = 2)] V 
^ i<0 ' 



Hence, 



dQ 
dP 



Y [E e -i u {Vi\X Tl = 1) + E e - i0J (Vi\X Tl = 2) 



The calculation of EQ-i w {Vi\XT 1 = k), (k = 1,2) is based on the branching process. Note that 
Vi = A(i + 1) + B(i + 1) + C(i + 1) + A(i) + B(i) = \U(i + 1)| + A(i) + B(i). Then 



X Tl = 1 



E u 



E W (|17(1)|+A(0)+B(0) P(l)JI Tl = l 



1 - a(0) - /3(0) 



1^(1)1 + 



l-a(0)-/?(0) 



1^(1)1 



X Tl =1 



1 - a(0) - /3(0) 
1 

1 - a(0) - /3(0) 



£ w (tf(i) x Tl = i 
(i,i,i),^ft/(i) 



x Tl = i 



Similarly, for i < — 1, 



EjVi 



X Tl = l 



E w [Vi 
1 



U(i + 1) 



l-a{i) -p(i) 
1 



1 - a{i) 
1 

1 - a(i) - /3(i) 



X Tl = l 
E u (U(i + l)\x Tl =l 
EjuQ.) X Tl = l) ■ N(0) • • • N(i + 1) 



(1,1,1),^ (U(l) X Tl =1 ■JV(0)---JV(i + l)). 



Hence, for % < — 1, 
X Tl = 1 



1 - a(0) - /3(0) 



(1,1,1), ^-^(f/(l)|x Tl =l) • AT(-i) • • • AT(1) 



By the same argument, 



And for i < — 1, 

X Tl =2 



X Tl =2 



1 



1 - a(0) - £(0) 



(1,1,1),^H7(1) X Tl =2 



1 



1 - a(0) - /3(0) 
By branching process, we have 

^(V(l) 



X Tl = 1 

(1,0,0) ^(£7(1) = (l,0,0)|X Tl = 1) 

+(0,1,0) -P U (U(1) = (0,l,0)|X Tl = 1) + (0,0,1) -P U (U(1) = (0,0,l)|X Tl = 1) 
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(1,0,0) 



Pi(0) 



Pi(0)+ 7 (0) 



1 - q(0) - 0(0) /VI- a(0) - (9(0) - 

7 (0) // pi(0)+ 7 (0) 



+ (0,1,0) 
7(0) 

Pi(0)+7(0)'pi(0)+ 7 (0) 



1 - a(0) - /3(0) / V 1 - a(0) - /3(0) 



+ (0,0,1) -0 



Pi(0) 



,0 



And, 



EjU{l) X Tl =2) =(0,0,1). 



Consequently, we get that 



§ = E(EMf=*)) 

i<0 V fc=l,2 7 

2 + /(l, 1, 1), E,>i f Efc=i, 2 ^ ( (i) ^ = fc) • JV(t) • • • N(l) 



1 - a(0) - /3(0) 
1 - q(0) - /3(0) 



□ 



4 The law of large numbers 

We are now ready to prove the law of large number with an explicit drift based on the branching 
structure by the method of "environment viewed from the particle" . 

Lemma 4.1 Under the law induced by Q <S> P u , the sequence {cJ(n)} is stationary and ergodic. 

Proof . Since Q is an invariant measure under the transition kernel K by Theorem ll.2| we obtain 
that the process {cJ(n)} is stationary under Q <%> P w . For the ergodicity, the proof is similar as 
Zeitouni ([H], page 207, Corollary 2.1.25), we omit the details. □ 

Proof of Theorem ] 1.31 The idea of the proof for the LLN is the same as Zeitouni (pj5], page 208, 
Theorem 2.1.9), however, we pay attention to the drift here. 

Define the local drift at site x in the environment uj, as d(x, u>) = E^(X\ — x). The ergodicity 
of {uJ(i)} under Q <g) P^ implies that: 

^ n— 1 n— 1 

— } d(Xk, oj) = — y d(0,uJ(k)) — > / d(0,ui) dQ, as n — > oo, Q <8> P^-a.s. 
n L — ' n I 

k=0 fc=0 
On the other hand, 

n n n 

X n = Y J {X i -X i _ l ) = ^(x i -Xi- 1 -d(X i - 1 ,u)}+Y i d(X i -. 1 ,u) 

i=l i=l i=l 
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Under P u , M n is a martingale, with \M n — M n _i| < 3 for uj £ 0. Hence, with ^ 
tr(Mi,--- ,M n ), 



and hence, iterating, E u (e^ Mn ) < e 3nA2 (this is a version of Azuma's inequality, see [?], Corollary 
2.4.7). Then Chebyshev's inequality implies that 



Mr, 



n 



0, P-a.s. 



Hence, 



Observe that 



lim — X n = I d(0,Lj) dQ = vp, 

n— >oo n 



Vp 



f d(0,u)dQ = Eq(X 1 ) 

n(a;)fpi(0) + 2p 2 (0)-p_i(0) 



E P 



Ep 



Q{n) 

Pi(0)+2p 2 (0)-p-i(0) 
l-a(0)-/3(0) 



E(T 1 \X 7i = l)+E(T 1 \X Tl =2) 
Moreover, by (|2.2|) and the proof of the density of Theorem \1.2\ we have 



E{Ti\X Tl = 1) = E P 
= E P 
= Ep 
= Ep 
and by the same argument, 



E w (l + ((2,2,1), £ 17(») 
l + ((2,2,l),^^([/(i) 



i<0 



X Tl = 1 
X Tl = 1 



i<0 



1 + ((2,2,1),^E W (?7(1) X Tl = lj • JV(0) • • • JV(i) 
pi(0) 7 (0) 



1 + ((2,2,1),^ 



i<0 



Pi(0)+ 7 (0)'pi(0)+ 7 (0) 



,0) -N{Q)-.-N(i)) 



E(Ii|X Tl = 2) = £ P 



1 + ((2,2,1), J] (0,0,l) • JV(0) • • • JV(i) 
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Therefore, 



E(Ti|X Tl = l)+E(Ti|X Tl =2) 
E P < Pl(0) 

E P 



2 + ((2,2,l),J] 



i<0 



7 (0) 



Pi (0) + 7(0)' Pi (0)+ 7(0) 



2+((2,2,l),^ 



7(0 



i>0 



px(i) + 7(i)' pi (i) +7(i) 



,i).jv(o)...j\r(t)) 

l) .N(i)---N(0)) 



where the last equation holds by the invariance of P. Consequently, the drift in this case can 
be written as 



Ep 



Vp 



Pi(0)+2p 2 (0)-p-i(0) 



l-a P (0)-^fo) UU; ( 2 + £)t>l U 1 ' 1 ' X )' ( pi(i)+7W PlW+7W 1 i ' N ® ' ' ' N ^ 



E P 



2 + U2, 2, 1), £. t >o ( ^g^y, ^^y, 1 ) • • • • iV(0) 



7p) 



□ 



Remark Consider the random walk in non-random environment defined in Proposition 11.11 
Then the drift vp reduces to E(X\) = pi + 2p2 — q. In fact, when the environment is not 
random, the drift in Theorem 11.31 can be written as 



vp = (pi + 2p 2 - q) 



l-a-P 



2+((l,l,l),(^,^,l)-En>l^ 



2 + /(2,2,l),f^,^,l).E n > 1 ^ 



By (pT22|h we obtain 



vp = (pi + 2p 2 



2 + 



(pi+7)(l-2Q-3/3) 



(1, 1, 1), (a(2 Pl + 3 7 ), P(2px + 3 7 ), 2/3 Pl + (1 - 2a) 7 



2 + 



(pi+7)(l-2a-3/8) 

2pi+37 

(pi + 2p 2 - q) x 



(2, 2, 1), (a(2p x + 3 7 ), /3(2px + 3 7 ), 2(3 Pl + (1 - 2a) 7 ) ) 



(pi+7)(l-2a-3/8) 



2pi+37 
(pi+7)(l-2Q-3/3) 



= pi + 2p 2 - g, 
as it should be. 



□ 



5 The general bounded jump case: (1, i?)-RWRE 

In this section, we consider (1, i?)-RWRE, in which the possible jumps to the right are 1, 2, • • • , R, 
where R is some fixed positive integer. In this case, the environment is an element 

w = {(g(w)„pi(w)*, • • • ,pr(uj) z ) :zeZ}e M(lf =-. n, 
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where A = { — 1, 1, ■ ■ ■ , R} is the set of all possible jumps. Following the idea in section [21 we 
can define a branching process with (1 + 2 + • • • + i?)-type species. In fact, for a jump from 

1 to i — 1, there are (1 + 2 + • • • + R) crossing-back ways, i.e., jumping from i — k\ to i for 
k\ £ {1, 2, • • • , R}, jumping from i — k<i to i + 1 for hi G {1, 2, • • • , i? — 1}, • • • , jumping from 

2 — /cr to i + i? — 1 for k R = 1 . 

First, we need the following lemma dealing with exit probabilities for this general case. 



Lemma 5.1 For n > 2 and 1 < j < R, we have 

(ei, [M(») + • • • + M(-n) • • • M(i)]( ej 



I*[(-n,i),i + j] 
with eR + i = and 



ej+i)> 



1 + (ei, [M(t) + • • • + M(-n) • • • M(i)]ei) 



/ Pi(i)H 



M(i) := 



1 



V 



Pfl-i(0+Pfl(0 Pit(Q \ 

<?(*) 9(0 





o / 



The proof is the same as Lemma 12.11 



□ 



Define for i < 0, 

U(i) = • • • , U R (i), U R+1 (i), • • • , U R+R -i(i), • • • , j7i +2 +...+K(i)l , 

where, for k = 1, 2, ■ ■ ■ , (1 + 2 + • • • + i?), C/fc(i) is the number of steps from i to i — 1 before time 
Ti with crossing-back from i — 1 to i, ■ ■ ■ , i — i? to i; i — 1 to i + 1, ■ ■ ■ , i — i? + 1 to i + 1; • • • ; 
i — 1 to i + -R — 1; respectively. For probabilities, let 

pm(i) = P^,[the walk jumping from i to i — 1, and finally jumping back 
from z — 1 to i before Ti], 

P(R)(i) = -fS[th e walk jumping from i to z 
from i — R to i before Ti], 
= T^[the walk jumping from i to 
from i — 1 to i + 1 before Ti], 

P(it+it-l) (?) = -P^fthe walk jumping from i to i — 1, and finally jumping back 
from i — R + l to i + 1 before T\], 

P( 1+ 2h \- R ) (i) = -P(!,[the walk jumping from i to i — 1, and finally jumping back 

from i — ltoi + ii — 1 before Ti]. 



— 1, and finally jumping back 
i — 1, and finally jumping back 
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Expressions of these probabilities can be calculated by using the exit probabilities as this is done 
in section [2j Firstly, we have 

P(l) + • • • + P(r) = -f^ [the walk jumping from i to i — 1, and finally 

jumping back to i before Ti], 
P(r+i) + • • • + P(r+r-i) = P^ [the walk jumping from i to i — 1, and finally 

jumping back to i + 1 before Ti], 

P((i +2+ ...+fl)-2) (i) + P((i+2+-+ J R)-i) (*) = Pi! [the walk jumping from i to % - 1, and finally 

jumping back to z + R — 2 before Ti], 

P(i+2h h-R) W = Pa, [the walk jumping from i to i — 1, and finally 

jumping back to i + i? — 1 before Ti]. 

Recall the exit probabilities, we obtain 

P(i) + ■ ■ ■ + P(R) = ?(*) • f£ -1 [(-oo,i - 1).*], 
P(«+i) + • • • + P(r+r-i) = g(i)--FS -1 [(-oo,i-l),i + l], 

P((i+2+...+fl)-2)(0 +P((i+2+...+fl)-i)(*) = ?(*) • Pi!" 1 ^-^^ - l),* + -R-2], (5.1) 

P(i +2+ ... + h)(«) = gCO-^K-oo.i-l^i + Jl-l]. 

Observe that 

P((l+2+...+fl)-2)(») : P((l+2+-+fl)-l)(») = Pfl-l(* - 1) : P(l+2+-+K)(« - !)• 

Thus by (15. If) . we get 



p {{1+2+ ... +R y 2) (i) = q(i) ■ P^K-oci- l),i + R-2] l{% l) -— -, 



P (( i +2 +... +J R)-i)W = ?(*) • Pi [(-oo,i - l),i + R-2] 



P(l+2+...+Jl)(*-l) 



Pfi-1 («-!)+ P(l+2+-+R) (* - 1) ' 
Using the same argument, we can get the expression of all the probabilities: 

p (1+2+ ... + fl)(i) = g(*)-i^ -1 [(-oo,*-l),< + 12-l], 
p (( i +2+ ... +H) _ 2) (i) = g^-P^K-oci-l),^^-^ PH-i^- 1 ) 

Pfl-l(* - 1) +P(l+2+-+R)0 - 1) 

p ((1+2+ ... +fl) _ 1) (i) = g(i)-Pr 1 [(-oo,i-l),i + i?-2] 



- 1) + P(i+2+-+R) (« - 1) ' 

p 2 (« - 1) 



p 2 (i - 1) + P(H+H_1 + 1) (» - 1) + • • • + p {i j +jR _i +J? _ 2 ) (* - 1) ' 
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= g(i)-i* -1 [(-oo,i-l),i + l] 

P(fl+fl-l+fl-2)(^ ~ 1) 

p 2 (* - 1) + P(H+fl-l+l) (* - 1) H H P(iJ +J R_l+iJ-2) (* - 1) ' 

= tfo-^K »,< u.,;- ^ h . / M ; i . n , 

P(2) (i) = gCtJ-i^-MC-oo.i-l),*] 



P(fl)(») = Q(i)-K [(-oo,»- l),t] 



pi(i 


- 1) + P(RH 




■■+P(HH 


-B-l) (* - 


1) 






P(h+i)(* - 


1) 






pi(z 


- 1) +P{m 


-!)(*-!) + • 


• • + P{m 


-n-i)(i - 


1) 






P(fl+Jl-1)(» 


-i) 







pi(z - 1) +P(H+i)(i - 1) H +P{R+R-i){i ~ 1)' 

Now we are ready to give the theorem about U(i). 

Theorem 5.1 Assume X n — > oo, P-o.s.. Then (U(i)) is an inhomogeneous multitype 

V J i<o 

branching process in R 1+2_l ^ R with immigration 

Pi(0) 



?7(1) = ei, ttratft probability 
U(l) = e 2 , with probability 

U(l) = en, with probability 
U(l) = eR+i, with probability 



U(l) = e2+3H KRi with probability 

U(l) = ei+2H \~r, with probability 



l-p (1) (0) p (Jl) (0)' 

g(g+i) (Q) 

1-P(i)(0) P( R) (0)' 

P(fl+fl-i)(0) 
1-P(i)(0) P( fl )(0)' 

P2(0) 

1-P(i)(0) p (Ji )(0)' 

P(l+2+... +J R)(0) 



1-P(i)(0) p ( fl)(0)' 

Pfl(0) 

1-P(i)(0) P( fl) (0)' 

and i/te following offspring distribution: 

P U (U(0) = («!,-•• ,«*,<),••■ ,0),17(1) = ei||l7(l)| = l) 
(ni H 1-%)! 



ui! • • • -urI 



: pS(o)---pS(o)pi(o), 



P w (l7(0) = («!,••■ ,«b,1,0,--- ,0),17(1) = e 2 ||C/(l)| = l) 

r) -pS(o)---p5(o)p ( h +1) (o), 



ui! • • • tiij! 
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/ (R+R-l)th 

P»[U(0) = (ui,--- ,u R ,0--- ,0, 1 ,(),••• ,0),C/(1) =e R \U(1)\ = 1 
(ui j h?i fl )! u 

P W (CT(0) = (ui,-- - ,« fl ,0,.-. ,0),C/(1) = e R+1 \\U{\)\ = l) 

= ( \ti--'-t7 )l p S (0) -" 1, c% (0)pa(0) ' 



P w (E/-(0) = («!,-■■ ,«*,(),••■ ,0,1),C/(1) =e 2+3+ ... +J R |C/(1)| = 1 

= (Ul u|!.'..^T )! p S (0) ' ' •«(°^ + - + K)(0), 
P w (f/(0) = (m, ■ ■ ■ , u R , 0, • • • , 0), U(l) = e 1+2+ ... +R \\U(l)\ = l) 

= ( " 1 ^:::t7 )l p s (0) -"^ (0)pa(0) - 

Furthermore, for i < 0, i/ie offspring mean matrix of the {—i + l)-th generation is: 

N(i) = (Nx(i) N 2 ), 



where 



P(i)(0 

l-P(i) (*) P(fl)(0 



P(i)W 

\ l-P(i) (*) P(K)(») 



P(fl)O) 

!-P(l) (*) P(K)(») 



P(fl)W 

!-P(i) W P(K)(») / 



(l+2+-+-R)xR 



and 



/ Zl t R-l 


Zl,R-2 


■ z 1A 




Ir-i 


Zr-1,R-2 ■ 


■ Zr-i,i 




Zl,R-l 


Zl,R-2 


■ z hl 




Zr~2,R-1 


Ir-2 


Zr-2,1 




Zl,R-l 


Zl,R-2 


■ z ltl 




Zl,R-l 


Zl,R-2 


• h 




V Zl,R-l 


Z\ y R-2 


■ z hl 


) 



No 



(l+2+-+K)x(l+2+-+-R-l) 

in which Z m ^ n is the zero matrix of dimension mxn, and I m is the identity matrix of dimension 
m x m. 



Remark Following the ideas in section [21 Theorem 1 5 . 1 1 can be proved analogously. We omit the 
details. 
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